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Abstract 

The diagrammatic linear response formalism for the Coulomb drag in 
two-layer systems is developed. This technique can be used to treat both 
elastic disorder and intralayer interaction effects. In the absence of intralayer 
electron-electron correlations we reproduce earlier results, obtained using the 
kinetic equation and the memory-function formalism. In addition we calcu- 
late weak-localization corrections to the drag coefficient and the Hall drag 
coefficient in a perpendicular magnetic held. As an example of the intralayer 
interaction effects we consider a situation where one (or both) layers are close 
to (but above) the superconducting transition temperature. Fluctuation cor- 
rections, analogous to the Aslamazov-Larkin corrections, to the drag coeffi- 
cient are calculated. Although the fluctuation corrections do not enhance the 
drag coefficient for normal-superconductor systems, a dramatic enhancement 
is found for superconductor-superconductor structures. 
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In a recent series of experiments the mutual friction between two parallel close electron 
systems [1-4], electron-hole systems [5] and normal metal-superconductor systems [6] was 
measured. In these experiments a current passing through one of the layers (the active layer) 
induces a current in the other layer (the passive layer), due to the frictional forces. If no 
current is allowed to flow in the passive layer, a potential bias is developed which cancels 
the frictional interlayer force. The trans-resistance (per square), or the drag coefficient, is 
defined as the ratio between the potential developed in the passive layer, to the current in 
the active layer, i.e. po = V pass i ve j '1 active- Measurements of pjy provide a direct measure of 
the frictional forces between the two layers. A few mechanisms can lead to frictional drag 
forces between the charge carriers of different layers. The carriers can interact via Coulomb 
forces [2,7], phonon exchange [3,8], or magnetic interactions of spontaneously created vortices 
(in the case of superconducting layers) [9]. The temperature dependence of the observed 
trans-resistance in electron [2] and electron-hole [5] systems [po oc T 2 ) was qualitatively 
explained by the Coulomb mechanism, although in order to explain the dependence of the 
trans-resistance on the distance between the layers a virtual phonon drag mechanism was 
suggested [3,8]. The present paper is limited to the study of the Coulomb drag mechanism 
only, however, our approach can be, to some extent, generalized to include other mechanisms 
of drag. 

Drag due to the Coulomb interaction between the carriers of spatially separated layers 
was considered more than a decade ago in Refs. [10,11]. The first direct experimental 
observation of drag and a comparison with theory was reported in Refs. [1,7]. The new 
experiments [2,5,6], as well as the exotic behavior of two layer systems in a strong magnetic 
field [12] have renewed recently the interest in the interlayer friction phenomena. 

The previous theoretical studies of the Coulomb drag were based upon the kinetic equa- 
tion [2,7,13] or upon the Mori memory-function formalism [14]. The latter formulation is 
sufficiently general to treat also disorder inside the layers. It was found that in dirty sys- 
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terns at a very low temperature, the temperature dependence of the drag coefficient changes 
from po oc T 2 to po oc T 2 lnT [f4]. In both cases, that of the kinetic equation and that 
of the memory-function formalism the approximations involved are uncontrollable, i.e. it is 
not clear what is the order of the corrections to the obtained results. It is also not clear to 
what extent these approaches can be used to treat quantum interference effects, the presence 
of an external magnetic field, or intralayer electron correlations. Therefore an alternative 
rigorous approach is needed. In this paper we employ Kubo linear response formalism in 
the framework of diagrammatic expansion in order to develop a controlled technique for 
calculating the drag in various setups. Besides a controlled way to obtain the old results, 
the diagrammatic approach provides a method to treat the intralayer interactions, quantum 
interference phenomena etc. We obtain a few new results, which can hardly be obtained (if 
at all) by the other techniques. 

The central quantity in the framework of the Kubo linear response is not the trans- 
resistance, pj), but rather the trans-conductance, (To- In order to obtain the experimentally 
measurable drag coefficient, one should invert the 2x2 conductance matrix. To lowest 
non-vanishing order in the interlayer interaction, one obtains 

where <T; = (e 2 /li)2Z);!/; is the conductance of the I th isolated layer (D is the diffusion 
constant, u-the density of states per spin). We emphasize the apparent analogy between the 
drag and the Hall coefficients since both of them are off-diagonal elements of the resistance 
matrix. As we shall see, this analogy goes very far and is extremely useful. The technical 
reason is that in both cases one has to calculate "three leg" vertices. The less formal reason 
is that the dragged current may be considered as a non-dissipative one since the applied 
voltage and the induced current are spatially separated. 

Before sketching the obtained results, let us list the limitations we have restricted our- 
selves to. We treat the drag effect only to the lowest non-vanishing order in the interlayer 
interactions, while intralayer interactions may be included to an infinite order. The interlayer 
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tunneling of electrons is not discussed here. We consider only thin quasi two-dimensional 
layers (the layer thickness should be less than the screening radius). Generalization to thick 
layers is not very painful, but requires some effort. We do not consider electron-hole sys- 
tems, but only two electron layers (possibly with different mobilities and Fermi energies). 
Only the spatially uniform d.c. drag effect is studied. The drag phenomenon may also occur 
between two types of carriers which are not spatially separated, e.g., in two band metals 
[15], heavy and light holes [16] etc. Although the diagrammatic approach is certainly useful 
for studying these examples, our results can not be applied directly. 

The previous theoretical studies of the drag phenomenon [7,8,13,14] emphasize the close 
relation between the drag coefficient and the imaginary part of the intralayer polarization 
operator (which is related to the dynamic structure factor). It was even proposed [17] 
that measurements of the drag coefficient may be a useful tool in studying the structure 
factor of strongly correlated systems. We show, that although the above mentioned relation 
works surprisingly well (e.g. it reproduces correctly the quantum interference effects and the 
response to an external magnetic field), it is not generally valid in the presence of intralayer 
electron-electron correlations. In general, as will be shown, the drag coefficient may be 
related to the non-linear susceptibilities of the layers. This suggests that the d.c. drag 
current is a result of rectification by the passive layer of an a.c. fluctuating electric field, 
created by the active one. In a particular case, when intralayer electron-electron correlations 
are absent, the non-linear susceptibility is reduced to the product of the diffusion constant 
and the imaginary part of the polarization operator, reproducing the earlier results. 

The role of intralayer correlations is demonstrated for the case in which one or both 
layers are close to the superconducting transition temperature. This part of our study was 
motivated by recent experiments on normal metal-superconductor systems [6]. We consider 
the contribution of spontaneously created fluctuating Cooper pairs to the drag. These lead 
to a divergent contribution (near the critical temperature) to the conductivity (Aslamazov- 
Larkin corrections [18]). We have calculated the contribution of these fluctuations to the 
non-linear susceptibility. It appears to be proportional to a divergent power of (T — T c ) (for 
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dirty two-dimensional films it is cx (T — T c )~ ). This divergence is not strong enough to 
enhance significantly the drag coefficient in normal metal-superconductor structures and, 
thus, fails to explain the experiments [6]. We found, however, that in the case of two identical 
metals, both close to the superconducting transition, the drag coefficient may be enhanced 
dramatically. The maximal enhancement (within the Ginzburg criterion for validity of the 
Landau-Ginzburg approach) is by a factor of <r - the dimensionless conductance in the 
normal state. 

The outline of the article is as follows. In Sec. II the diagrammatic technique for the drag 
trans-conductance is developed. The main result of this section is Eq. (6), which establishes 
the relation between the trans-conductance and the non-linear susceptibility. In Sec. Ill 
the non-linear susceptibility, T, is calculated for (clean and dirty) normal metals (without 
intralayer electron-electron interactions). In this case T ~ DSII (II is the polarization 
operator), and the old results are reproduced. The dependence of the trans-resistance on the 
mobilities, the interlayer distance and the temperature is listed in Sec. Ill C and summarized 
in Table I for completeness. In addition to the known results [14], we point out a novel 
linear temperature dependence of pjj in a certain range of temperatures, which takes place 
in very dirty systems. In Sec. IV we study the quantum interference (weak-localization) 
corrections to the trans-conductance. Due to a cancellation of a certain class of diagrams, 
the expression derived in the memory-function formalism [14] appears to be applicable 
(with the proper renormalization of the diffusion constant). In Sec. V we show that in the 
absence of intralayer interactions the Hall drag coefficient vanishes (again in agreement with 
the memory-function formalism). Sec. VI is devoted to the study of fluctuation corrections 
to the drag coefficient near the superconducting transition. We conclude with a discussion of 
the results and a brief summary. The appendix is devoted to the calculation of the screened 
electron-electron interactions in two layer systems. We do not compare the theory with the 
existing experiments, mainly because it can be found in the literature (see eg. Ref. [4]). The 
new results obtained here are waiting for an experimental confirmation. 



II. KUBO FORMULATION OF TRANS-CONDUCTANCE 



Repeating the standard derivation of the Kubo formula [19] for the case in which an 
external field is applied to layer 1 and the induced current is measured in layer 2, one 
obtains for the drag trans-conductance 

= [ ' dte iat (\ji\Q,t),Ji(Q,0)\). (2) 

Here i, j = i, y\ Q, fi are the wave vector and the frequency of the external field, S is the area 
of the sample and J{ is the i th component of the current operator in the I th layer. Diagrams 
corresponding to Eq. (2) include two separate electron loops with a vector (current) vertex 
on each one of them, coupled only by the interlayer interaction lines. To second order in the 
interlayer interaction there are two different diagrams, depicted in Fig. 1. (The first order 
diagram vanishes when the external field wave vector Q is zero. That is the case we are inter- 
ested in.) We shall focus hereafter on the d.c. trans-conductance = <Jp(Q = 0, SI — v 0). 




FIG. 1. Two diagrams, contributing to the trans-conductance to second order in the interlayer 
interaction. The full lines with arrows are the electron Green functions. The wavy lines represent 
interlayer interactions. Full dots are vector (current) vertices. 

Analytically, the two leading order diagrams (Fig. 1) may be written in a symmetric 
form 

■ T 

Od = Q ^ r i(£ w " + ^m^n)^i(q^n^n 4- Q m )U(q,u n + Q m )U(q,u n ), (3) 
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where U(q,co) is the interlayer screened Coulomb interaction (see Appendix A) and the 
vector r i ( 2 )(^*, ^>i, UJ2) is the "three leg" object given by the two diagrams depicted in Fig. 2 
(the factor 1/2 in Eq. (3) is included in order to prevent double counting). 




FIG. 2. Diagrams defining the "three leg" vertex, r(g, ^i,^) — r( — q, — —tjJ\). 

In Eq. (3) the usual Matsubara technique with fi m — 27rimT is employed [20]. After 
summing over the boson frequencies, one should perform an analytical continuation 
to a real value of H, and finally the limit Q — > should be taken. The sum over io n is 
done by a contour integration in the complex to plane along the contours shown in Fig. 3a. 
Since the integrals over the small circles cancel the contributions from points Lo n — and 
Lo n — —Q m [21], one finds that 



D 2 



E J- l dz coth 7^ V M * + On, z)Yi(l z, z + tt m )U(l z + Q m )U(l z) 

q 

= -^E /l^ coth ^ [r^-^^-n^rr+^-n,^) (4) 
-rj+-(£ w + n ra ,w)r 2 J " + (q,w,a,' + n m )] t/ + (<f,w + n m )tr (<?», 

where the symbols + and — indicate from which part of the complex plane the analytical 
continuation is performed. The way of analytical continuation of r + ~(^*, a^,^) is illustrated 
in Fig. 3(b). Note that only the middle part of the contour (cf. Fig. 3(a)) contributes to the 
Eq. (4). The upper and lower parts of the contour do not contribute, because in the limit 
0^0 
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f^(q,LO,oo) oc ^(SH - |^ n± (<?» = 0. 



(5) 



where H^^q^io) is the retarded (advanced) polarization operator (the two terms in the 
middle part of the last equation correspond to the two diagrams of Fig. 2). As a result, 
r ±± (g, lo + 0,tj) oc 0, and do not contribute to the d. c. trans-conductance. 





CO 










U> 2 





(a) 



(b) 



FIG. 3. (a) The contour of integration in the complex u> plane, employed to perform the sum 
over the Matsubara boson frequencies in Eq. (3). (b) The way in which the analytical continuation 
of r(<j, £Ji, ^2) — » T + ~ (q,bj,uj) is performed. 

We perform now the analytical continuation, Q m — » 0, in Eq. (4) and take the limit 
SI — > 0. The result is 



i fO 

Jd = 167rTS^J-c 



-r;+- (q, lo, Lo)T 3 2 - + (q, lo, lo) \U+(q, lo) 



(6) 



Here we have used the fact that U (q,u) — (ll + (q,u;)j . Applying Onsager relations to 
<Tj)(H)) one can easily show that 



T+-(q,UJ,UJ,H) = T-+(q, 



(7) 



where H is an external magnetic field. 

Eq. (6) provides the relation between the d. c. trans-conductance, cr^J, and the vertices, 
ff^(q,iv,uj). The latter are the non-linear susceptibilities (namely, rectification) of the 
electron gases in an external a.c. field. The physical interpretation of this fact is the 



following. If the passive layer is exposed to an external space and time dependent scalar 
potential, <f)(q,oj), then a d.c. current, 

J 2 = f^(<f,^M<?>)| 2 , (8) 

is induced, as a result of the second order rectification effect. In the drag problem the 
random field, <f)(q,io), originates from the thermal fluctuations of the electron gas in the 
active layer. In the presence of the external driving field, E\, the correlator of the non- 
equilibrium component of the random field is 

<^)~(£-iY-)^|tf + | 2 , (9) 

where rae(faj) is the Bose distribution function. Substitution of Eq. (9) in Eq. (8) leads to 
Eq. (6). 

The next section is devoted to the calculation of the non-linear susceptibility, 
r + ~ lo, to) in the absence of the intralayer electron-electron interactions and an exter- 
nal magnetic field. It will be shown that, to 0(1 /(kpl)), T in this case is given by 

f + -(q,LO,Lo) = f- + (q,LO,Lo) = e— $ffl + (<f,fai), (10) 

where D = D(Q — 0,0 — 0) is the d.c. diffusion constant of the electron gas. 

In this particular situation, employing Eqs. (1) and (6), one obtains for the d.c. trans- 
resistance the result of Refs. [7,13,14]: 

p% - ~ ^£<?V T ^^3n^>)3n+( 9 >)|u+( 9 >)| 2 . (n) 

We stress, however, that Eq. (10) and hence Eq. (11) are proved below (up to order 
0(l/(kjpl))) only for the case in which electron-electron correlations within each layer are 
absent. For interacting electrons (e.g. for the drag between two superconductors) one 
should use the general expression, Eq. (6). The formal reason is that the electron-electron 
interactions corrections cannot be reduced, in general, to a renormalization of the vertices 
in the diagrams for f, Fig. 2 (see also Fig. 5 bellow). 
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III. DRAG BETWEEN NORMAL METALS 



We turn now to the calculation of the non-linear susceptibility, r + ~ (<f, fa;, fa;), for normal 
metals, without intralayer electron-electron interactions. It will be shown that in this case 
Eq. (10) may be obtained within controlled approximations. 

Recall that r (q, lo\ , lo 2 ) is given by the two diagrams depicted in Fig. 2 and the analytical 
continuation should be performed according to Fig. 3(b). The expression corresponding to 
Fig. 2 has the form 

f(q,w u w 2 ) = T£Tr + G^G r _^M_ n -^Y (12) 

where I is the current operator, Q tn is the Green function and the trace is taken over the 
exact quantum states of the disordered system. Next we perform the sum over the fermionic 
energies, e„, and the analytical continuation of T according to Fig. 3(b). This is done in a 
similar way to the one employed for the boson frequency sum in Eq. (3) (cf. also Fig. 3(a)). 
Then e.g. for the first term on the r.h.s. of Eq. (12) one has 

— / 'del tanh — - tanh Tr (c+C" J IQ+ \ (13) 

4vri J \ IT IT ) \ ' E+ ^ 2 J V ; 

+ h I dt ( tanh ^w- - tanh 2?) Tr {^^ + Jst + ^ } , 

where denotes retarded (advanced) single electron Green functions. There are no con- 
tributions from e n > — fa, 2 > and e n < < since the trace of three Green functions of 
the same type vanishes. Performing the analytical continuation according to Fig. 3(b), one 
obtains 

+ {q,u -> -q,-w} . 

Employing in advance the fact that Tr [- - ■] does not depend explicitly on e (for epr ^> 1), 
and recalling that in the absence of a magnetic field Q~ — (<? £ + )*, one obtains 

f +- w , fa>) = ^STr \$; F G: F+ jGt F+ J\ + {?-, " - -I • (15) 
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Special care should be taken in the presence of an external magnetic field (see our discussion 
of the Hall drag coefficient in Sec. V). Further evaluation of the non-linear susceptibility, 
r +_ (<f, a;, a>), should be performed separately for different values of the momentum and 
energy exchange, q,co. 



A. The diffusive case 

To calculate r + ~ (<f, to, to) (see Eq. (15)) in the diffusive regime, i.e., when to <C 1/r and 
q <C 1// (t and / — u^r are elastic mean free time and path correspondingly), one should 
employ the impurity diagrammatic technique (see e.g. Ref. [20]). We use here its simplest 
version with short range impurities only. The extension to the general case is straightforward 
and we shall comment on it at the end of this section. The dominant diagram is depicted 
in Fig. 4. 



FIG. 4. The leading order (in (M)" 1 ) dia gram for T + (q,w,w) in the diffusive regime. The 
analytical properties of the Green functions are indicated by + (retarded) and — (advanced) signs. 
Dashed ladder denotes a Diffuson [20]. 

The corresponding analytical expression is given by 



2e, 



7rr Dq 2 — iiv . 



ii 



In the last line we took an advantage of the fact that lot\ ql <C 1 and expanded up to the first 
non-vanishing order in these parameters. The momentum sum in the last expression should 

be performed with some care. The naive evaluation results in: ^ . . . = 27T Svt 2 Dq and the 

p 

two terms on the r.h.s of Eq. (16) cancel each other. The cancellation occurs because the 
asymmetry between electrons and holes is ignored by the linearization of the spectrum near 
the Fermi energy . Without such asymmetry the electrons' drag and the holes' drag compen- 
sate each other completely. A very similar situation arises in the calculation of the Hall coeffi- 
cient of dirty metals in the framework of the diagrammatic technique [22,23] . In order to take 
into account the asymmetry we approximate 1 y^ j v p (—v p q) ... — —Svq I de p (ej? -\- e p )/m . . . 

J — CO 

p 

and obtain ^ ... — 2tt Svr 2 Dq (l -\- -^—^j ■ As the leading order does not contribute to T + ~ , 

p 

one should keep track of terms of the order (e^r) -1 . (For a non-parabolic dispersion relation 
€f should be replaced by (de p /dp) 2 (2d 2 e p /dp 2 ) 1 at the Fermi energy. ) Substituting the 
result for ^ . . . into Eq. (16), one obtains 



I (q,uj,uj) — e 



9? ^ 
isv- 



(Dq 2 ) 2 +w 2 J ' 

which has precisely the form of Eq. (10), since in the diffusive regime [24] 



(17) 



Da 2 

n + (g» = 2Si/ - - q — . (18) 
Dq 2 — ilo 

We conjecture that Eq. (10) is valid if diagrams of the form of Fig. 5(a) are considered. 
Here (and in the next section) this conjecture is proved for specific examples. The diagrams 
of the form of Fig. 5(b), which cannot be reduced to a renormalization of the vertices, make 
Eq. (10) invalid. However, to leading order in (cpr) -1 , these diagrams do not contribute. It 
will be shown in the next section that they do not contribute to weak-localization corrections 
as well. We conclude thus that Eqs. (10) and (11) are valid up to 0(-^) in the diffusive 
case. One can easily convince himself that the inclusion of long range scattering will not 
change this conclusion. Indeed, the small angle scattering renormalizes the current vertex 
and one of the scalar vertices, leading to a redefinition of the diffusion constant only. 
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(a) 



(b) 



FIG. 5. The general structure of diagrams for the non-linear susceptibility, T + ~(q, to, to). The 
black areas represent a vertex renormalization by impurity and intralayer interaction lines. The 
shaded triangle is an irreducible "three leg" vertex. 

In order to perform the calculation of the drag coefficient one should employ the screened 
interlayer interaction in the diffusive regime which is given by (cf. Appendix A) 

1 7re 2 q D\q 2 — ico D 2 q 2 — ito 



U + (q,Lo)- 



(19) 



S KiK^sinhqd D\q 2 D 2 q 2 ' 

where the divergences at small q are cut off at D^{2)1 2 % lJ /( fi: i(2)^)i here ^1(2) are the 
Thomas-Fermi momenta. 



B. The Ballistic case 

Next we consider the ballistic regime where either q > 1//, or lo > 1/r. In this case 
the dominant contribution to r +_ comes from diagrams without any internal impurity lines 
(each Green function is still impurity decorated). Employing the approximation 

GT+MQSp) = 2r30 e - ^(p) »i 27tt<5(c + lo - e p ), (20) 

substituting it in Eq. (15) and performing the momentum integration, one obtains for 
lo; q 2 /2m <C t F 
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f(q,u,u) = e?^\2Sv—0(v F q-u)\, (21) 
e F 1 v F q 1 

where 0(x) is the Heaviside step function. As in the diffusive case the two terms on the r.h.s 
of Eq. (15) cancel each other in the leading order. To obtain Eq. (21) the non-linearity of 
the dispersion relation should be taken into account. Note, that Eq. (21) has the form of 
Eq. (10) again. Indeed, for the ballistic system (at small q and lo) [19] 

3Tl + (<?» = 2Su— 8(v F q -lo). (22) 
v F q 

The screened interlayer interaction in the ballistic case is given by (cf. Appendix A) 

U + {*M = \ ™\, (23) 
b K1K2 sinned 

After establishing the validity of Eq. (10) (and therefore Eq. (11)) for the diffusive and 
the ballistic cases we shall list results for the drag coefficient for various temperatures and 
mobilities. 



C. Non— interacting drag trans-resistance for different temperatures and mobilities 

In order to understand how the drag coefficient depends on the temperature, the mobil- 
ities and the distance between the layers it is instructive to look at the relevant parts of the 
(<?,w) plane (see Fig. 6). 

We distinguish between two parts of the plane 

Ballistic for lo < vj?q a > 1//; 

(24) 

Diffusive for lo < 1/r q < 1/1. 

The frequency integration according to Eq. (11) is cut off exponentially at lo % T, while the 
momentum summation is restricted (also exponentially) by q % 1/d due to the interlayer 
interaction potential. We thus distinguish between ballistic, / ^> d, and diffusive, / <C d, 
samples. All the experiments [1,2,5] until the present day were done on ballistic systems. 
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FIG. 6. Domains in the (q,w) plane. Dashed area represents (schematically) the regions where 

3n+(<2» ^ 0. 

For ballistic systems (/ ^> d) the dominant contribution to the trans-resistance, Eq. (11), 
comes from the ballistic part of the (q,co) plane (apart from the exponentially small range 
of temperatures [14]). Substituting Eqs. (22) and (23) into Eq. (11), one obtains [2,7,13,14] 

% 7T 2 ((3) T 2 1 1 

e 2 16 e F1 e F2 (n\d){K,2d) (k F1 d)(k F2 d) 
For diffusive systems (/ <C d) the entire contribution comes from the diffusive part 
of the (q,io) plane. Substituting Eqs. (18) and (19) into Eq. (H), one obtains for 
T <C min{r _1 ,To} (where To = D min k 2 } jd arises due to the divergence of the in- 
teraction potential at small momenta, cf. Eq. (19)) 

= A£Ll!_i T ° 1 1 (26) 

PD en2e F1 e F2 n 2T(K 1 d)(K 2 d)(k F1 l 1 )(k F2 l 2 ) ■ ( ' 

For larger temperatures, T ,T ^> t -1 , the energy integration is dominated by the region 
lj % r _1 , the result is 

~ % IT 1 1 (27) 

e 2 £F\£F2 (/Ci<i)(/C2<i) (kFih)(k F 2h) 

(in order to determine the exact numerical prefactor one should know the behavior of 

3II + (q,tij) for lo fa t -1 .) Actually, in this case there is an additional contribution due 
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to two-dimensional plasma modes [25]. In fact, the contribution of the plasmons in this 
region of temperatures appears to be dominant. 

For extremely dirty samples, / -C d/Viai -C d, a situation where r _1 > T ^> T may 
occur. The drag coefficient in this case is 

= s_xf(3) _rn_ i i _ (28) 

e 2 8 e Fl e F2 {^id)(K 2 d) (k^l^k^h) 
In order to obtain Eq. (28) the full expression for the screened interlayer interaction, 
Eq. (A2), should be employed. 

The dependence of the drag coefficient on the temperature and the interlayer distance is 
summarized in Table I. 



TABLE I. Temperature and interlayer distance dependence of the drag coefficient, /j£j, for 
different mobilities. 





T < min{T ,r- 1 } 




T > mincer" 1 } 


1 < d 


d- 2 T 2 \nT 




d- 2 T a 


I > d 




d -4 T 2 b 





a For T > r 1 the contribution of the plasma modes is dominant [25]. 

b For exponentially small temperatures T < r _1 exp{-3(//d) 2 /4f (3)}, p D ~ d- 2 T 2 \ogT [14]. 



If) 



IV. WEAK-LOCALIZATION CORRECTIONS TO THE DRAG COEFFICIENT 



In this section we calculate the weak-localization corrections to the non-linear suscep- 
tibility, r +- (^*, lj : u), and to the drag coefficient, pp. We restrict ourselves to the diffusion 
case, T -C t _1 and / <C d, as this is the case where the localization corrections can be 
essential. The drag coefficient, without the localization corrections is given by Eq. (26). 
Three types of diagram, depicted in Fig. 7(a)-(c), may contribute to the weak-localization 
corrections to r + ~ lj : co). 




(a) 



(b) 



(c) 




(d) 



FIG. 7. (a)-(c) The weak-localization corrections to the non-linear susceptibility, r + (q,u:,u:). 
The zigzag lines are Cooperons [26], the shaded rectangle (d) is the Hikami box [26,27]. 

The first two diagrams, Fig. 7(a), (b), belong to the family represented in Fig. 5(a). As we 
discussed above, such diagrams, lead to the result given by Eq. (10). Indeed a straightforward 
calculation of diagrams Fig. 7(a), (b) yield 



2D(0,0)« C 



■2Si. 



(i)( g » 9 2 ) 2 + u 2 



(29) 



where the renormalized diffusion coefficient D(q,Lo) for a two-dimensional system is given 
by [26] 
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D(q,u>) = D (l--^- In— ). (30) 

V 7rfci?/ LOT J 

Eq. (29) was verified up to the first order in SD = D(q,co) — D; one should understand 
Z)(0,0) as Z^OjT^ 1 ), where t^ is a phase relaxation time [23,24], 

The diagram in Fig. 7(c), is of the type of Fig. 5(b). It can not be reduced to the 
product of the d.c. diffusion constant and the polarization operator, like Eq. (10). It is 
easy to check, however, that due to Hikami-like cancellation [26,27] this diagram does not 
contain divergent factors like lnr/r^ or Incor. It is not surprising then (this can be checked 
by a direct calculation 1 ) that the diagram in Fig. 7(c) is smaller than the value given by 
Eq. (29) by a factor of (kpl)^ 1 . We neglect, thus, the contribution of this diagram. As 
a result, Eq. (29) accounts correctly for the first order weak-localization corrections to the 
nonlinear susceptibility. This expression has again the form of Eq. (10), ensuring the validity 
of Eq. (11). Employing Eq. (30), we obtain for the weak-localization corrections to the drag 
coefficient 

^ = -_i-]nJ---i-]nJ- ) (31) 
pD 7rk F ih 2Tti 7rk F2 h %Tr 2 

where pjj is given by Eq. (26). The last equation should be compared with the weak- 
localization corrections to the conductivity of a two dimensional system [23,26] 

S<7 \ ^ Ta . 

— = - -T-M-- 32 

The main difference is that in the case of drag the logarithmic singularity is cut off by T 
instead of As usual, one expects to have a dependence on a weak magnetic field [23]. 

In the absence of spin-orbit scattering we predict a positive drag magneto-trans-resistance. 
The characteristic scale of the magnetic field is given by H c — T/(e_D) (in the case of 
conventional weak-localization it is H c — l/(e_Dr^) [23], which is typically much smaller). 



■technically this calculation is very similar to the calculation of the weak-localization corrections 
to the Hall coefficient [22,23], where two Cooperon diagrams are also canceled. 
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The intralayer interaction corrections to the magneto-resistance have typically the same 
characteristic scale of magnetic fields [24], making the observability of the weak-localization 
corrections to the drag coefficient problematic. 



V. HALL DRAG 

Next we consider the two-layer system in a perpendicular magnetic field. In order to 
calculate a X pj we follow the usual treatment of the Hall conductivity in disordered systems 
[22,23]. The external electric field (in the active layer) and magnetic field (in both layers) 
are introduced by E — iVlA-i and H — i |& X yl 2 J respectively. The Hall current is propor- 
tional to E x H = (A 2 (k ■ Ai) - k(A 1 ■ A 2 )). It is convenient to choose a gauge where 
(k ■ A 2 ) — 0. The first order corrections (oc H) to the non-linear susceptibility, T + ~ (<f, oj, oj), 
are given by the three diagrams depicted in Fig. 8. In order to calculate the contributions of 
these diagrams we expand the Green functions and retain only the contributions which are 
proportional to A 2 (k ■ A±). Performing the momentum summation, one should keep track 
of terms up to order l/(kjrl), since the leading order terms are canceled out. The procedure 
is similar to the one that was employed in the derivation of Eq. (17). 




FIG. 8. The three digrams contributing to the first order corrections in the magnetic field to 
the non-linear susceptibility, r +_ (q, w, to). 
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For example, for the diagram Fig. 8(a) one obtains: 

L a ~ ^v p (eA 2 ■ v p _ k/2 )Q+ +UJ (p)g- (p - q- k/2)Q~ +UJ (p- k)Q~ +U) (p) (33) 

p 

« 2^^(eA 2 ■ v p - k/2 )(q- v p ){k-v p )Gt(p){G;(p)f -+ 2e(A 2 ■ q)k\su2TTT 5 v F (l - ~rz ~ yj )■ 
? 8 i(t F T} 

In the same manner the other diagrams result in ~ e(A 2 ■ q)k-^Sf27^r 5 Vf(—2 — 2 ( e ^ T )2 )> 

L c — 0. Some care should be taken while performing the above calculations. The point is 

that in the presence of a magnetic field the contributions from Q~Q~Q + and Q + Q~Q + (see 

Eq. (14)) are not complex conjugated and should be calculated separately. (Note that each 

of these terms separately is not gauge invariant and depends on the particular choice of 

momentum parameterization (gauge). Only the sum of the two, which is independent of the 

momentum parameterization is gauge invariant.) Adding all the terms together, one obtains 

to first order (in the magnetic field) corrections to the non-linear susceptibility 

ST+-(q,L0,L0) = e— \qx S] —W+tq^). (34) 
tp L 1 m 

Substituting this in Eq. (6) and employing Onsager relations, Eq. (7), we obtain 

a*J = (w c n + ^r 2 )^, (35) 

where co c — eH/m is the cyclotron frequency. 

In order to calculate the Hall drag coefficient one has to invert the 4x4 conductance 
matrix (with the components lx,ly,2x,2y), keeping only the leading order in both the 
interlayer interaction and the magnetic field. The result is 

Pl v - 0. (36) 

The same result can be obtained by direct application of Eq. (11) to the xy component of the 
trans-resistance [28]. Thus, it has been proven that the off-diagonal part of Eq. (11) is valid 
as well as the diagonal one (in the absence of the intralayer electron-electron interactions). 
The fact that ppf — is closely related with the cancellation of the diagram Fig. 8(c) 
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(L c — 0). The physical interpretation of this observation is that the magnetic field does 
not influence the interlayer interaction but rather affects the "free motion" (diffusion) of the 
electrons. As a result the Hall drag current is given by a sum of two obvious contributions: 
the Hall component of the dragged current in the passive layer and the current dragged by 
the Hall component in the active layer (cf. Eq. (35)). This statement, however, may be not 
correct for higher orders of interaction corrections. 



VI. DRAG NEAR THE SUPERCONDUCTOR TRANSITION 

In this section we consider the Coulomb drag when one (or both) layer is close to (but still 
above) the transition to the superconducting state. Namely, we shall consider the vicinity of 
the transition temperature, where the fluctuations of the superconducting order parameter 
become important. These fluctuations lead to the Aslamazov-Larkin (AL) [18] corrections 
to the conductivity. For dirty (1/r <C T c ) quasi two-dimensional films the AL corrections 
have a remarkably simple form [18] 

Sa AL ^^—-, £ = ^1_^ < 1 5 (37) 

where T c is the transition temperature and a is the phenomenological constant in the 
Ginzburg-Landau equation (in the framework of BCS theory a — tt/8). 

Our aim is to find the fluctuation corrections to the drag coefficient. We restrict ourselves 
to dirty superconductors, T c <C 1/t. The extension to the clean case is straightforward, 
but the fluctuation corrections are less pronounced there. As the intralayer (attractive) 
interactions play a crucial role, one cannot use Eq. (11) and has to employ the general 
expression, Eq. (6). (Application of Eq. (11) with the AL corrections to the Sn + leads to 
a wrong result.) The most divergent (as T — > T c ) corrections to the vertex, T, are given by 
the diagram depicted in Fig. 9(a). 
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9^2 



(a) 




(b) 



FIG. 9. (a) The fluctuation corrections to T(q, ^2)1 near the superconducting transition. 
The bold zigzag lines represent propagators of fluctuating Cooper pairs [18]. (b) The triangular 



This diagram may be considered as the non-linear susceptibility of the Ginzburg-Landau 
order parameter in an external a.c. field (in a complete analogy with the electron non-linear 
susceptibility given by the diagram, Fig. 4). The corresponding analytical expression is 



8V(q,io U Lo 2 ) = 8T K $ ~ £ Xm)B Xm , Xm+UJ2 K(k, Xm + io 2 )C t K(k, Xm 4- io x )B Xm+UJuXrr , 



(38) 



where the factor 8 is a result of the charge of a Cooper pair and the spin summation. In the 
above equation K(k,Xm) is a propagator of the fluctuating Cooper pairs, 



C t = eSv—2Dk 

k rp 



(39) 



is the vector triangular vertex [18] (see Fig. 9(a)) and B XUX2 is the scalar triangular vertex, 
Fig. 9(b). Integrating over a contour with three cut lines, \ — 0j ~ ^lj ~ ^2, in the complex 
X plane and performing analytical continuation according to Fig. 3(b), we obtain 

6f+-(q,co,Lo) = 8 £ ^ jH d X (™th^ - coth^) K~(k, X )K + (k, x ) (40) 

", :-.. xJ +{q,u->-q,-u;}. 



\K + (k -q, x + ^)B^B^ X - K~(k -q, x + J)BtZ^B. 
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The terms with K + K + K + and K K K cancel the contributions of the unphysical poles at 
X — and x — ~ LO ou 4- m the close vicinity of the transition, £ < 1, one may approximate 



cotn — coin- 



2T c lo 



(41) 



2T 2T x(x + u)' 

If we perform now the integration in Eq. (40), employing the standard expressions for the 
fluctuation propag ators, K*, [18], we shall obtain an exact cancellation, i.e., r + — (after 
the term with q, lo — > — q, —lo is taken into account). This is again due to linearization of 
the spectrum (see Sec. Ill A). The same problem arises in the calculation of the fluctuation 
corrections to the Hall coefficient [29-31]. It was shown [31], based only on the gauge 
invariance of the Ginzburg-Landau equation, that the fluctuation propagator, which takes 
into account the electron-hole asymmetry, has the form 



Sv 



a , , (a i d In T c 
-Dk Tlx l- T - — 



+ e 



(42) 



where (i is a chemical potential. Note, that — —edK~ x jdk. The form of the propagator, 
Eq. (42), suggests that the scalar vertex, B cx dK~ x j dx, should also contain the asymmetry 
factor, dlnT c f d/i, [32]. An explicit calculation of the triangular scalar vertex, Fig. 9(b), 
leads to the following results 



-Hxi 



X2) 



dlnT c 



T c Dqi-i( Xl -X2)' 2 dp 
«■ -Hxi + Xi) _^i_d\nT c 



(43) 



Next we substitute Eqs. (39),(41)-(43) into Eq. (40) and perform the energy integration. 
Leaving only the most divergent (as £ — > 0) terms, we obtain to first order in the asymmetry 
factor, dlnT c /d[i, 

Sa 2 dlnT c LoDq 2 ^ 1 T + + T_ 



6T + ~(q,u,u) = e2Dq- 



d[i (Dq 2 ) 2 + w 2 "? T + T_ (T + + T_) 2 + (aLo/T c ) 2 ' 



(44) 



T ± ^-£(*±|) 2 + £ . 
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We shall see that the dominant contribution to the drag trans-conductance comes from 
an extremely small interlayer momentum exchange, Dq 2 fa (T — T c )/(kcI) <C (T — T c ). The 
intralayer momentum exchange, k, has a characteristic value given by Dk 2 fa lo fa (T — T c ). 
Since q <C k the summation over k in Eq. (44) may be easily performed, leading (for the 
two dimensional case) to 



ST + (q, a;, lo) — e- 



2Dq 



3 In Tc \ 

an ) 



2>V 



oDq 2 



(Dq 2 ) 2 +L0 2 



■(—, 

\2T c eJ <t«£ 2 



(45) 



^) = -i^ln(l + x 2 ) 

and <7 S — 2Dv is the dimensionless conductance of a film in a normal state. Eq. (45) (cf. 
with Eq. (17)) solves the problem of finding the fluctuation corrections to the non-linear 
susceptibility (for Dq <C ui) in dirty superconductors. In the next subsections we shall 
incorporate it with the general formula, Eq. (6), in order to find the drag trans-conductance 
between a normal metal and a superconductor and between two identical superconductors. 
Let us stress that Eq. (45) can not be obtained by a substitution of the fluctuation corrections 
to D(q,co) into Eq. (16). The latter procedure leads to a wrong conclusion, ST + ~ cx e _1 , 
instead of £r +_ cx e~ 2 . 

The validity of the perturbation theory is restricted by the condition 8<Jal *C <7 s , i.e. 
e ^> Gi, where Gi — a/(27r<r s ) is the Ginzburg parameter for dirty two-dimensional super- 
conductors [30]. In the temperature range of 

(46) 



Gi < £ < VGi 



the fluctuation corrections to the non-linear susceptibility become dominant, i.e. ST ^> T 
(see Ref. [30] for a discussion of the analogous situation in the case of the Hall coefficient). 



A. Superconductor— Normal Metal Drag 



As we already mentioned above, the fluctuation corrections to the drag coefficient come 
from extremely small energy and momentum interlayer exchange, lo fa (T — T c ) <C 1/r; 
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Dq 2 % (T — T c )/(/cef) <C 1/t. Because of this one should employ expressions derived for 
the diffusive regime, Sec. Ill A, irrespective of the ratio between d and /. Substituting 
Eqs. (17), (45) and (19) into Eq. (6), we obtain the fluctuation corrections to the drag 
trans-conductance 

where the subscripts s and n refer to the superconductor and the normal metal respec- 
tively; k = min{>c n ; k s }. Comparison with Eqs (26) and (37) shows that up to logarithmical 
multipliers 

^£!~^~_L (48 ) 

Since the corrections to the drag coefficient are given by 

(49) 



they are at most of the same order of magnitude as pjj in the normal state. Note that, due 
to the integration over the full (q,io) phase space, Sp s £~ n < pj) even if ST ^> T. Hence, we 
think that the fluctuation corrections to the Coulomb drag coefficient cannot explain the 
recent experiments [6] (for an alternative explanation see [9]). 

B. Superconductor— Superconductor Drag 

Next we consider the drag effect between two identical superconductors. We confine 
ourself to the temperature region Gi <C e <C y/Gi, where ST ^> T. Substituting Eqs. (45) 
and (19) into Eq. (6) and performing the integrations, we obtain 

e 2 l-ln2 T? ln(«.d) a 

In the present case the divergence appears to be extremely pronounced, oc £~ 3 . As a result 
the fluctuation corrections to the drag coefficient may exceed essentially it's normal state 
value. For e fa Gi the enhancement is Sp"^' / pn ~ <t s ^> 1. The normal state conductance, 
<t s , cannot be arbitrarily large since the dirty film condition, T c < l/i, should be fulfilled. 
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VII. DISCUSSION OF THE RESULTS 



In this section we summarize and discuss the main results of this paper. First, the regular 
diagrammatic treatment of the interlayer Coulomb drag to the lowest non-vanishing order 
in the interlayer interaction has been developed. The method that was developed may be 
useful in other situations where two kinds of mutually interacting carriers are important 
(e.g. two band materials). Non-Coulomb mechanisms of drag may be easily treated within 
the same scheme (to second order in the interlayer interaction). 

It has been shown that the calculation of the drag trans-conductance may be reduced 
to the calculation of the non-linear susceptibility (rectification) of the electron gas, Eq. (6). 
This fact reflects the physical mechanism of the drag phenomenon. Namely, the d.c. drag 
current is a result of the rectification by the passive layer of the fluctuating electric field, 
created by the active one. In special cases, the non-linear susceptibility may be reduced 
to the product of the diffusion constant and the polarization operator, Eq. (10), which 
reproduces the earlier result, Eq. (II), derived using the kinetic equation [7,13], and the 
memory function formalism [14]. We conjecture that this reduction is valid, when the 
correlations between the electrons inside each layer may be neglected. The above conjecture 
was checked by a direct calculation of the non-linear susceptibility for diffusive and ballistic 
systems. It proved to be valid even for the off-diagonal trans-conductance (Hall drag) and 
the weak-localization corrections. In the last two cases the validity of Eq. (10) is not obvious 
a priori and reveals itself through a subtle cancellation of certain classes of diagrams (Figs. 
7(c) and 8(c)). For completeness the dependence of pd on the temperature, the mobilities 
and the interlayer distance in the absence of intralayer interactions are listed in Sec. Ill C 
and summarized in Table I. As a byproduct we point out the novel linear temperature 
dependence of the drag coefficient, taking place in extremely dirty samples, I < \Jdj k. 

Next a system where the intralayer electron-electron correlations are of crucial impor- 
tance is studied. As an example, metals just above the transition to the superconducting 
state were investigated. Fluctuations of the superconductor order parameter cause cor- 
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rections to the non-linear susceptibility. The dragged current may be enhanced due to 
spontaneously created unstable Cooper pairs. The non-linear susceptibility is found to di- 
verge as (T — T c )~ 2 close to the superconductor phase transition. Note that Eq. (10) which 
works perfectly well for the non-interacting case, gives a wrong result, cx (T — T c ) _1 . Having 
the result for the susceptibility, one can easily calculate the fluctuation corrections to the 
drag trans-conductance between a normal metal and a superconductor and between two su- 
perconductors. In the first case the correction to the drag trans-conductance appears to be 
proportional to (T — T c ) _1 . As a result the drag trans-resistance (Eq. (49)) is not enhanced 
significantly. Apparently the fluctuation corrections can not explain the recent experiments 
[6] (for an alternative explanation see Ref. [9]). In the case of two superconductors, however, 
the divergence of the trans-conductance appears to be much stronger, cx (T — T c )~ 3 . As a 
result the trans-resistance may be enhanced by a factor of <r 3 ^> 1 (<r s is the normal state 
dimensionless conductance of the film). We believe that the last prediction may be checked 
experimentally. The example of superconductors shows that the relation between the drag 
coefficient and the polarization operators, Eq. (f 1), may be incorrect for strongly correlated 
systems. 
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APPENDIX A: SCREENED INTERLAYER INTERACTION 



In this Appendix we treat the screened Coulomb interaction in two layer systems. To 
this end we employ the standard random phase approximation (RPA) [19]. 




FIG. 10. Screened Coulomb interactions in the RPA. The dashed and full wavy lines are the 
bare and the screened interactions correspondingly. 

The Dyson equation (see Fig. 10) in the RPA is most conveniently written in a matrix 
form 



14 U 

u v 2 



1 ' v° u ^ 



u° v° 



I 



v° u° 
u° v° 



lit 
IF 



14 U 

u v 2 



(Al) 



Here V° — 2?re 2 / (Sq) and U° — V°e~ qd are the bare intralayer and interlayer interactions 
correspondingly; Vi and U are the screened intralayer and interlayer interactions, Tli(q, uj) 
is polarization operator of the I th layer. Since the tunneling between the layers is neglected 
the off diagonal elements of the n matrix are equal to zero. Solving the RPA equation and 
employing the bare interaction value, we obtain for the screened interlayer interaction 

,2?re 2 



2n+(q,w)n+(5,w)sinh(^)- 



Sq 



Sq 
27re 2 



+ n+(«, w ) + n+(«, w )), 



(A2) 



In order to obtain Eqs. (19) and (23) we substitute Jl + (q,Lo) — ISvDq 2 j (Dq 2 — ito) for 
the diffusive system, or Il + (q,Lo) = 25V for the ballistic one and retain the leading order in 
> 1 (k = Aire 2 i/ is the inverse Thomas-Fermi screening radius). 
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